
function [r] = R_Solve(q)
global  D E  ;

% Solve Er=-Dq
r = -E\(D*q);
